{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### 1. Uploading data\n",
    "The code requires command ``akdensity`` to be installed. \n",
    "Type ``ssc install akdensity`` if it is not installed yet."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "use \"https://github.com/worldbank/Stata-IE-Visual-Library/raw/develop/Library/Density%20plots/Shaded%20k-density%20functions/data.dta\", clear"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### 2. Create variables for the graph"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "sort beta_\n",
    "gen rank = _n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### 3. Calculate parameters for graph"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "qui egen p98=pctile(beta_), p(98)\n",
    "\n",
    "qui sum beta_, det\n",
    "local mean = round(`r(mean)', .00001)\n",
    "local median = round(`r(p50)', .00001)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### 4. Creating the graph"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "akdensity0 beta_, gen(x) at(beta_) bwidth(.0005) //akdensity0 comes from the user-written package \"akdensity\"\n",
    "\n",
    "qui sum beta_, d"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "twoway area x beta_ if rank>15 & beta_<(`r(p10)'), color(gs14) || ///   \n",
    "\tarea x beta_ if beta_>`r(p90)' & rank<980, color(gs14) || ///      \n",
    "\tarea x beta_ if rank>15 & beta_<(`r(p5)'), color(gs9) || ///    \n",
    "\tarea x beta_ if beta_>`r(p95)' & rank<980, color(gs9) || ///    \n",
    "\tline x beta_ if rank>15  & rank<980, lcolor(black) || ///    \n",
    "\t(pcarrowi -20 .00299 310 .00299, lcolor(cranberry) lpattern(dash) msize(zero)) || ///    \n",
    "\t(pcarrowi -20 `r(mean)' 310 `r(mean)', lcolor(gs7) lpattern(dash) msize(zero)) || ///    \n",
    "\t(pcarrowi -20 `r(p50)' 310 `r(p50)', lcolor(gs7) lpattern(dash) msize(zero)), /// \n",
    "\tlegend(off) ///\n",
    "\txtitle(\"2SLS Coefficient from baseline model\" \" \") ///\n",
    "\tytitle(\"Density\" \" \") ///\n",
    "\txlabel(0 \"0\" .00299 \"NQ=.00299\" `r(p50)' \"Median=`median'\" `r(mean)' \"Mean=`mean'\" .015 \".02\", angle(45)) ///\n",
    "\tylabel(none) ///\n",
    "\tbgcolor(white) graphregion(color(white))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### 5. Exporting the graph"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "> (gs7) lpattern(dash) msize(zero)),  legend(off)  xtitle(\"2SLS Coefficient from baseline model\" \" \")  ytitle(\"Density\" \" \")  xlabel(0 \"0\" .00299 \"NQ=.00299\" `r(p50)' \"Median=`median'\" `r(mean)' \"Mean=`mean'\" .015 \".02\", angle(45))  ylabel(none)  bgcolor(\n",
      "> white) graphregion(color(white))\n",
      "\n",
      ". noi gr export C:/Users/Saori/.stata_kernel_cache/graph${stata_kernel_graph_counter}.png, width(600) replace        \n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAlgAAAG0CAIAAABhas7vAAADAFBMVEUAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAACzMPSIAABS+klEQVR42mP4PwpGwSgYBaNgFIxgwPB/FIyCUTAKRsEoGMFgtCIcBaNgFIyCUTCiwWhFOApGwSgYBaNgRIPRinAUjIJRMApGwYgGoxXhKBgFo2AUjIIRDUYrwlEwCkbBKBgFIxqMVoSjYBSMglEwCkY0GK0IR8EoGAWjYBSMaDBaEY6CUTAKRsEoGNFgtCIcBaNgFIyCUTCiwWhFOApGwSgYBaNgRIPRinAUjIJRMApGwYgGoxXhKBgFo2AUjIIRDUYrwlEwCkbBKBgFIxoM0orwyJEj/0fBKBgFo2AUjALag0FaESoqKv4fBaNgFIyCUTAKaA9GK8JRMApGwSgYBSMajFaEo2AUjIJRMApGNBitCEfBKBgFo2AUjGgwWhGOglEwCkbBKBjRYLQiHAWjYBSMglEwosFoRTgKRsEoGAWjYESD0YpwFIyCUTAKRsGIBqMVIdXAt1sP/w8+8P3bm/8kgvfvfv4fBaNgFIyCEQNGK0LqgL8/fx1kM/1+/+n/wQSePXt08mD9fxLB5N6LT5+++D8KRsEoGAUjA4xWhNQBfz59Pchq8n7/6f+DCVy6dKmwsPA/iaCgoODy5cv/R8EoGAWjYGSA0YqQOmC0IhwFo2AUjIIhCkYrQuqA0YpwFIyCUTAKhigYrQipA0YrwlEwCkbBKBiiYLQipA4YrQhHwSgYBaNgiILRipA6YLQiHAWjYBSMgiEKRitC6oDRinAUjIJRMAqGKBitCKkD/v78dULJ5/PZ6/8HE7h9+3ZjY+N/EkFjY+Pt27f/j4JRMApGwcgAoxXhKBgFo2AUjIIRDUYrwlEwCkbBKBgFIxqMVoSjYBSMglEwCkY0GK0IR8EoGAWjYBSMaDBaEVIH/P3564Jjyrc7j/8PJvD48ePJkyf/JxFMmjTpyZMn/0fBKBgFo2BkgNGKkDpgdPvEKBgFo2AUDFEwWhFSB4xWhKNgFIyCUTBEwWhFSB0wWhGOglEwCkbBEAWjFSF1wGhFOApGwSgYBUMUjFaE1AGjFeEoGAWjYBQMUTBaEVIHjFaEo2AUjIJRMETBaEVIHTBaEY6CUTAKRsEQBaMVIXXA35+/joo7fTx+6f9gAtevX6+urv5PIqiqqrpx48b/UTAKRsEoGBlgtCIcBaNgFIyCUTCiwWhFOApGwSgYBaNgRIPRinAUjIJRMApGwYgGoxXhKBgFo2AUjIIRDUYrQuqAvz9/XfbJ+37/6f/BBJ4/fz5z5sz/JIKZM2c+f/78/ygYBaNgFIwMMFoRUgeMbp8YBaNgFIyCIQpGK0LqgNGKcBSMglEwCoYoGK0IqQNGK8JRMApGwSgYomC0IqQOGK0IR8EoGAWjYIiC0YqQOmC0IhwFo2AUjIIhCkYrQuqA0YpwFIyCUTAKhigYrQipA0YrwlEwCkbBKBiiYLQipA74+/X7IW7LDwfP/h9M4MqVKyUlJf9JBMXFxVevXv0/CkbBKBgFIwOMVoSjYBSMglEwCkY0GK0IR8EoGAWjYBSMaDBaEY6CUTAKRsEoGNFgtCIcBaNgFIyCUTCiwWhFSB3w9+evaxHlg+3Q7ZcvXy5YsOA/iWD+/PkvX778PwpGwSgYBSMDMPwflGDIVYSj2ydGwSgYBaNgiILRipA6YLQiHAWjYBSMgiEKRitC6oDRinAUjIJRMAqGKBitCKkDRivCUTAKRsEoGKJgtCKkDhitCEfBKBgFo2CIgtGKkDpgtCIcBaNgFIyCIQpGK0LqgNGKcBSMglEwCoYoGK0IqQNGK8JRMApGwSgYomC0IqQO+Pfn77+fv/79+ft/MIG/f//++fPnP4ngz58/f/8OLo+MglEwCkYB7cBoRTgKRsEoGAWjYESD0YpwFIyCUTAKRsGIBqMV4SgYBaNgFIyCEQ1GK0LqgL8/f91IaRyEh24vX778P4lg2bJlL0cP3R4Fo2AUjBgwWhFSB4yuGh0Fo2AUjIIhCkYrQuqA0YpwFIyCUTAKhigYrQipA0YrwlEwCkbBKBiiYLQipA4YrQhHwSgYBaNgiILRipA6YLQiHAWjYBSMgiEKRitC6oDRinAUjIJRMAqGKBitCKkDRivCUTAKRsEoGKJgtCKkDhitCEfBKBgFo2CIgtGKkDrg789fP1++/fPl2//BBH7+/Pn58+f/JIJPnz79/Pnz/ygYBaNgFIwMMFoRjoJRMApGwSgY0WC0IhwFo2AUjIJRMKLBaEU4CkbBKBgFo2BEg9GKkDrg789ft/O7BuGh22vXrv1PIli7du3L0UO3R8EoGAUjBoxWhNQBo6tGR8EoGAWjYIiC0YqQOmC0IhwFo2AUjIIhCkYrQuqA0YpwFIyCUTAKhigYrQipA0YrwlEwCkbBKBiiYLQipA4YrQhHwSgYBaNgiILRipA6YLQiHAWjYBSMgiEKRitC6oDRinAUjIJRMAqGKBitCKkDRivCUTAKRsEoGKJgtCKkDvj789f3u49/v/v4fzCBr1+/vnnz5j+J4PXr11+/fv0/CkbBKBgFIwOMVoSjYBSMglEwCkY0GK0IR8EoGAWjYBSMaDBaEY6CUTAKRsEoGNFgtCKkDvj789eDhhmD8NDt7du3/ycRbNu27eXoodujYBSMghEDRitC6oDRVaOjYBSMglEwRMFoRUgdMFoRjoJRMApGwRAFoxUhdcBoRTgKRsEoGAVDFIxWhNQBoxXhKBgFo2AUDFEwWhFSB4xWhKNgFIyCUTBEwWhFSB0wWhGOglEwCkbBEAWjFSF1AHUrQk9PTw8Pj/8Ug9GKcBSMglEwCgiC0YqQOoBaFeHXr18ZkMDMmTP/UwBGK8JRMApGwSggCEYrQuqAvz9/fbl46+cLkk+4RgMMDAwMDAynTp06fvw4IyMjAwNDV1fXf3LBp0+fnj4leY//06dPP3369H8UjIJRMApGBhitCAcRgNR858+f/////////799+wYR+T8KRsEoGAWjgGZgkBayI7AiLCgoYGBgMDc3/48Ejh07xsDAwMjI+H8UjIJRMApGAW3AaEU4WAAjIyMnJ+d/DBAUFMTAwJCUlPR/FIyCUTAKRgENwGhFSB3w9+evRz0LyT50W1hYmIGBYcmSJf+xAUZyB0hfvny5b9++/ySCvXv3vhw9dHsUjIJRMGIAOcUrHcCQqwgpXDXKyMjIxMT0HweYOHEiAwODl5fXfxLB6KrRUTAKRsEoIAhGK0LqAEoqwp07dzIwMGzZsuU/bsDGxsZI+kzhaEU4CkbBKBgFBMFoRUgdQElFyMXFxcrK+h8vOHv2LAMDAy8v739SwGhFOApGwSgYBQTBaEVIHUBJRcjAwNDY2PifEJCUlGQgcaZwtCIcBaNgFIwCgoC0gpVuYORUhC4uLgwMDC9evPhPCFy9epWBgUFNTe0/0WC0IhwFo2AUjAKCYLQipA4guyJkY2MTFhb+TxwQEREhaaZwtCIcBaNgFIwCgmC0IqQOILsiZGBgmDp16n/iwP79+xkYGDw9Pf8TB0YrwlEwCkbBKCAIRitC6gDyKsJ169YxMDDcvn37P9GAkZQ9haMV4SgYBaNgFBAExBapdAZDriL8+/PXx6MXfjwiPNWHDDg5OfFsH8QKli9fzsDAsGDBgv9EgPfv39+/f/8/ieDevXvv37//PwpGwSgYBSMDMPwflGDIVYTkAQYGBuInCOGAhYWF1OpzFIyCUTAKRgEuMFoRDhj49u0bAwPDjBkz/pMI3r17x8DAIC4u/n8UjIJRMApGAcVgtCIcMGBoaMjAwHDv3r3/pANRUVEGBoYvX778HwWjYBSMglFAGaBKRfjk3q1Pv/5TEwy5ivDvz19Pp68m6dBtNjY2BqKXvaCBb9++sbCw6Onp/ccLXr58eeTIkf8kgiNHjrwcPXR7FIyCUTBiAJkFMSrYV8DAwC5qGNi69ASVasQhVxGSsWqUi4uL+I0QmGD+/PkMDAylpaX/cYPRVaOjYBSMglFAEFClIvz/69OjCztm57pJ81KpRhwJFSEDA8PkyZP/UwCYmJiY8K6aGa0IR8EoGAWjgCCgTkUIB78+PbqwbkI8okZcd+UlORXisK8Ia2pqGBgY7ty5858CMGnSJAYGhokTJ/7HAUYrwlEwCkbBKCAIqFwRfn93ff/sykBDUXYGdlHDwEA3aV4G2bAVN0mtDId9Rejj40PwxgliADMzMyPuQ9dGK8JRMApGwSggCKhTEULrP21eBgYGXmm33Nn7r7/7/v////////+6Od2avfTgf9LAsK8ImZiYrKys/lMMVq5cycDAEBkZ+R8bGK0IR8EoGAWjgCCgSkUIWSxjHz9h3YVHGFODX/f2524ldQRw2FeEjIyMV69e/U8NwMTExIBj9eloRTgKRsEoGAUEAfYClESAbfvEr0+vnsN6haSD4V0R3r17l4GBYdeuXf+pAXbv3s3AwCAgIPAfA4xWhKNgFIyCUUAQMPynAthXwFCw7z8q+LUrC1OQaDC8K8Ls7GwGBoaPHz/+pxLw9vZmYGCIi4v7jwpGK8JRMApGwSggCCiqCG8vClRXV1dXl+Zl4JVWRwXSvAy6E8ke+xtyFeHfn7/e7TpO5IZ6FxcXXl7e/1QFjNhupXj9+vWNGzf+kwiuX7/++vXr/6NgFIyCUTAyAHrRSRr4dXF5eXl5ebAug25wOSqYsPTMo1//yQVDriIkCQgJCRkbG/+nKli9ejUDAwMbG9v/UTAKRsEoGAWkAMoqQii4srp89ZX/1ATDuyJkZGT8+/fvf2qDnJwcBgaGgoKC/6NgFIyCUTAKiAYM/8kHT87vPP/k/////5+c34kVQGTJAcO4Irxx4wYDA8O+ffv+0wCoq6szYAyQjoJRMApGwSjAAygpNOFrZPYVMGAFEFlywJCrCP/+/PVi8RZi5ghzc3MZqLpSBhk8e/aMAemGppcvX54+ffo/ieDUqVMvRw/dHgWjYBSMGEBJRUhDMOQqQuJXjaqrqzPhPSCUQlBXV8fAwHDixIn/o6tGR8EoGAWjgAhAg4rw+9MLByk8c3s4V4S8vLwSEhL/aQk4ODhYWFj+j1aEo2AUjIJRQASgTkX469GuSreWo/////91skGWgYGBgYFdufLgx//kgmFcEbKysoaGhv6nJVi6dCkDA0N7e/toRTgKRsEoGAUEAVUqwo/bktllk7c9gbBEs3a9/fVye5ao55In/8kEw7UifPXqFQMDw9atW//TGEA6haMV4SgYBaNgFBAEVKkID5ayx27++P///1/7CtjZC/b9+v////99BSNpsQyRFeGGDRsYGBiePCG7hUAs2LdvHwMDg729/WhFOApGwSgYBfgBVSrCfQUMWbt+/f//61i1KHvyto///////3FzLOmXTsDBcK0IAwIC6LbnXUhISFlZebQiHAWjYBSMAvyAKhXhizVhvIalS9c3OrJD68G357od2QNWvPhPJhiuFaGOjo6GhsZ/egEVFZX8/PwlS5b8JwWMVoSjYBSMghEFqFIRQhbLSEtLu1Xue/L///////cV8LrPu/nrP7lguFaEXFxcaWlp/+kFjh8/np+fz8DAsGrVqv9Eg9GKcBSMglEwogB1KkKqgyFXEf79+ev1ur0EN9QzMDDU19f/pxd4+fLlgQMHWFlZGRgYeIk+5vvixYsvRzfUj4JRMApGDKBWRfj96YWD62e3lKMA8g8gHXIVITHg0KFDDDQ7XA0/gNSFdJueHAWjYBSMgiEEqFIR/rrYrcuAeROTesvR/2SCYVkRZmZmMjIyPn369P9AACEhIQYGBi4urv+jYBSMglEwCpAAVSrCkw2iuu3nvv+nHhiWFaGQkBAvL+9AVYT////n5eVlYGBwcnL6PwpGwSgYBaMABqhSER4sZSd/yyBWMOQqwr8/f71evw//HCEbG5u0tDQ9K8KXL19eunTpPxKwtLRkYGBwcHD4jxuMzhGOglEwCkYUYPhPBfBxW7Js1q63/6kHhlxFSMyqUXl5+ZqaGnpWhJewnSzDzMzMgPeqptFVo6NgFIyCEQXwFYhEgyurE6z42UUNA7PKUcAIWixDsCKEHK62aNGiAa8Inz9/zoB34cxoRTgKRsEoGFGAKhXh0RZ1rGAELZYhWBHOmTOHgYHh5MmTA14R/v//f+LEiQwMDHv37v2PDYxWhKNgFIyCEQWoUhFSHwy/ijAmJoaNje3hw4eDoSL8//8/MzMzIyPjf2xgtCIcBaNgFIwoQL2K8PvTCweXToAOhz65duXlr//kg+FXEbKysnJycj548GCQVIRVVVUMDAzBwcH/McBoRTgKRsEoGFGAShXhnTVhsgy80tK80BsnDpayyxbse/ufXDD8KkJGRkZra+vB0yP8//8/Ozs7I7ZO4WhFOApGwSgYUYAqFeHHbcm8jrNv/vq/D3710ttdWaKO8+7/JxMMv4qQgYGhtbX18ePHg6ciPHr0KAMDg6qq6n9UMFoRjoJRMApGFKBKRXiwlD1r16///5EqQlQ2yWCYVYRnzpxhYGDYtGnToKoI////LyAgwICxlWK0IhwFo2AUjCiAXgiSBQ7CLuZFqvxG2H2Ef3/+erFoC64N9d3d3QwMDLdu3aJzRfjy5ctTp079xwsYGBi4ubn/I4FTp069HN1QPwpGwSgYMYAqFeHHbcm8jtOvfodXhN/vrUuWhV5NSBYYchUhfmBkZMTPz//48ePH9K0IiQGurq4MDAx37tz5PwpGwSgYBSMSUKUi/P//7b4C6GIZXQs3dWleBgbZsDUUFK3DrCJ0d3c3MDB48uTJIKwI////z8nJKScn938UjIJRMApGJKBSRfj////vTy/sgN7DNGHpiXuUHcE9zCpCQUHB6OjoZ8+eDc6KsKioiIGBobOz8/8oGAWjYBSMPEC9ipCqYMhVhH9//nq3+zjWOcKfP38yMDDMnDmT/hXh69evb9y48Z8IICAgICIi8v///////9+4ceP169f/R8EoGAWjYGQAyivCXy+v7JjdEu+mrq6urq5uFpjVMnv/rU+//lMEhlxFiGfVKGSlzMGDB+lfERJcNQoHixcvZmBguHHjxv/RVaOjYBSMghEGKKsIfz3anKnMjnrcdrybNC8Du3IpBdvph1dFmJSUxMDA8Pjx48FcEf7//19QUBCyfHS0IhwFo2AUjChAUUV4Z7Y1u2HDoZfo3b9fL093O/IGrHjyn1wwnCpCd3d3SUnJx4O+Iuzq6mJgYLhz585oRTgKRsEoGFGAkorw5nQz9e6L/7GCX+fa1a1nk71udDhVhDIyMpqamk+ePBnkFeH///85OTn19PRGK8JRMApGwYgClFSEsF2D2AF+WQJgOFWEPDw89fX1T58+HfwVoba2NgMDQ8FoRTgKRsEoGElgtCKkDsBVEb5584aBgWHZsmXPnj17NugrwmvXrjEwMBSMVoSjYBSMgpEERitC6gBcFeGiRYsYGBjOnTv37NmzZ4O+Ivz//7+EhMRoRTgKRsEoGFGAwopQN7gcFwjWHa0I//+vqqoSFBR89OjRs2fPng2FinDatGkFBQV79uz5PwpGwSgYBSMDUFIRHm1Rxw9ajv4nEwy5ivDvz19Pp6/C3FDPx8dnaGg4UBXhy5cvjxw58p8U8OPHD319fW1t7f+jYBSMglEwMgAlFSENwZCrCHEBKSkpGxubp0+fPnv27BndK0LyADs7OwPG3UyjYBSMglEwXMEgLe+GTUXIwMCQkZHxDAaGREW4YcMGBgaGZ8+e/R8Fo2AUjIIRAEYrQhqCZ8+eMTAwtLW1PYOBIVER/v//n42NzdfX9/8oGAWjYBSMADBaEVIH/P356+Oxiz8ev/iPBBYuXMjAwLBo0aJnMEDnivD9+/f379//TyK4f/++uLg4ExPT/1EwCkbBKBgBYLQipA7Aumq0o6ODk5Nz/fr1z2CAzhUhGatG////X1BQMGPGDAYGhnfv3v0fBaNgFIyC4Q5GK0LqAKwVYWZmpoaGxoYNG57BwFCpCM+ePcvAwGBoaPh/FIyCUTAKhjsYrQipA7BWhCoqKiYmJkOxIrx8+bKRkREzM/P/UTAKRsEoGO5gtCKkDsBaEbKxsYWFhQ3RijA7O5uBgeHJkyf/R8EoGAWjYFiD0YqQOgCzInzx4gUDA0Nubu4QrQi/f//OyMgYGxv7fxSMglEwCoY1GK0IqQMwK8Lp06czMDC0trYO0Yrw////RkZGLCws/0fBKBgFo2BYg9GKkDoAsyKcMmUKPz//hAkTNgzZirCtrY2BgeHPnz//R8EoGAWjYPiC0YqQOgCzItTT01NTU5s0adLQrQgho7sVFRX/R8EoGAWjYPiC0YqQOuDvz1+PuhciH7qtrKysp6c3sBXhy5cv9+7d+59EsHfv3pcvX/7////////KysqjR8yMglEwCoY3GK0IaQUYGRmDg4MHtiKkHEyfPp2Hh+fv37//R8EoGAWjYJiC0YqQJuDt27cMDAxZWVmTJ08e0hXhy5cvGRgYVqxY8X8UjIJRMAqGKRitCGkCjh8/zsDA0NLSMtQrwv///3t4eCgpKf0fBaNgFIyCYQpGK0LqgL8/f325dPvny7f/////////vHnzODk5Z86cObAV4adPn56RfpvS06dPP3369B8G5s2bx8DA8PXr1/+jYBSMglEwHMFoRUgdgLZqNC0tTUFBYc6cOQNbEVK4ahQCPn/+zMDAICws/H8UjIJRMAqGIxitCKkD0CpCOTk5c3Pz+fPnD4OK8P///1xcXKysrP9HwSgYBaNgOILRipA6AK0ilJeX9/PzW7x48eRhURFu27aNgYHhypUr/0fBKBgFo2DYgdGKkDoAuSL88eMHAwNDamrqsmXLhkdF+P//f05OTk1Nzf+jYBSMglEw7MBoRUgdgFwRHj58mIGBob6+fjhVhM3NzYyMjP9HwSgYBaNg2IHRipA6ALkizM/PZ2RknDFjxnCqCG/evMnAwLBy5cr/o2AUjIJRMLzAaEVIHYBcEXp7e3NwcMyePXs4VYT///+XlJTk4OD4PwpGwSgYBcMLjFaE1AHIFaGrq6uiouL8+fOHWUU4YcIEBgaG9+/f/x8Fo2AUjIJhBEYrQuqAvz9/3a+bBjl0W15e3tvbe9GiRQNeEb58+XLbtm3/SQRbt259CTt0Gxn8/PmTmZk5IiLi/ygYBaNgFAwjMFoRUhn8/PmTgYEhJSVl2bJlywa6IqQ68Pb25ubm/j8KRsEoGAXDCIxWhFQGR44cYWBgqKmpWbZs2bJhVxHu37+fgYHh8OHD/0fBKBgFo2C4gNGKkMogOTmZgYEBsmR02bCrCP///8/Hxze6oXAUjIJRMJzAaEVIHfD356/v9578fveprKxMXFwcsmR02UBXhN++fXv7FnoOOPHgzZs33759+48DQJbM/Pjx4/8oGAWjYBQMCzBaEVIHwFeN+vj4mJubL1iwYNmyZcsGuiKk7qpRCPj8+TMjI6Opqen/UTAKRsEoGBZgtCKkDoBXhOLi4v7+/osXL162bNmy4VgR/v//39ramoFhkKacUTAKRsEoIBUM0uJsiFaETzYfYGBgyMrKWgYDw7IifPz4MQMDQ25u7v9RMApGwSgY+mC0IqQOgFSE+/pmMTAwtLW1LYOBYVkR/v//39HRkYmJ6f8oGAWjYBQMfTBaEVIHQCrC+XnVAgICs2bNWgYDw7UivHbtGgMDQ0hIyP9RMApGwSgY4mC0IqQOgFSENT4RWlpac+fOXQYDw7Ui/P//v6WlpYSExP9RMApGwSgY4mC0IqQOgFSEdrySurq6kMPVIGAYV4SXL19mYGBITk7+PwpGwSgYBUMZjFaE1AGQitCIkTc+Pn4ZEhjGFeH////5+fkZRpePjoJRMAqGOBikpdiQqwj//vx1KLRAkpG9trZ2GRIY2Irw5cuXa9eu/U8iWLNmzUtsh25jgt27dzMwMFRXV/8fBaNgFIyCIQtGK0KqgYkTJ7KwsMAPV4OAga0I6QCsra1ZWFj+j4JRMApGwZAFoxUh1YCTk5O4uDjySpllI6AiPHfuHAMDAxn9zlEwCkbBKBgkYLQipBpQVla2s7ODH64GAcO+Ivz//39kZKSdnd3/UTAKRsEoGJpgtCKkDvj785cMN39saPjSpUuXIYGBrQh//vz5+fPn/ySCz58///z58z/R4OTJkwwMDJs3b/4/CkbBKBgFQxCMVoTUATfPXzrIatIZmboMFQxsRUjrVaNwoKamJiQk9H8UjIJRMAqGIBitCKkDVi9acpDVZFFh7TJUMEIqwgsXLjAwMNja2v4fBaNgFIyCoQZGK0LqgPrS8oOsJstKG5ehghFSEf7//9/b25uHh+fPnz//R8EoGAWjYEiB0YqQOiDIw/sgq8mK8uZlqGDkVITv379nZmY2MjL6PwpGwSgYBUMKjFaE1AGK4pIHWU1WV7UtQwUjpyL8////lClTGBgYuru7/4+CUTAKRsHQAaMVIRXA8+fPuRiYD7KarK3pXIYKRlRF+P//f29vbwYGBj8/v/+jYBSMglEwRMBoRUgFsG7dutGKEA5sbW0ZGBg4OTlfvHjxfxSMglEwCgY9GK0IqQBaW1shQ6OjFSEEdHd3MzAwMDIyZmZm/h8Fo2AUjILBDUYrQiqAkJAQcyPjLfYxKydOX4YKBrYifPny5bJly/6TCJYtW/aSuEO38YP6+noGBgYGBgZubu4VK1b8HwWjYBSMgkEJRitCKgApKSk/P7/FixcvwwADWxEOBtDV1SUoKMjAwMDAwMDLy3vixIn/o2AUjIJRMJjAaEVIKXj27BkDA0NOTs4ybGC0IoQDd3d3BgYGBgYGJiYma2trY2NjaWlpBgYGDg4ONjY2BgYGNjY2RUXF169f/x8Fo2AUjAI6gtGKkFKwYsUKBgaGyZMnL8MGJo9WhKhg165doaGhioqKampqampqWlpatra2U6ZMeffunbi4OAMDAwMDAwcHh4yMzJo1a/6PglEwCkYB7cFoRUgpyM7OVlFRmTtz1vIFi5YtXrIMFQxsRfjnz5/fv3//JxH8+vXrz8AdELNkyRIxMTEGBgYGBgYWFhYmJqawsLDZs2dPnjz59u3b/0fBKBgFo4DaYLQipBSoqam5ubktmTVndNUodcH169fnzp2rqKjIzMzMAANMTEwsLCzv37//PwpGwSgYBVQCoxUhReDFixcMDAxZWVnL5y0YrQhpB96/f//kyZPjx49nZGQwMTExMTEJCwv/HwWjYBSMAmqA0YqQIrB69WoGBoZJkyaNVoR0A3/+/FFQUGBgYJCRkSkpKfk/CkbBKBgFlIHRipAikJeXp6SkNGfOnNGKkM7gxo0b6enpDAwM4uLi69ev/z8KRsEoGAXkgtGKkCKgqanp4uKycOHC0YpwQMCFCxf09PQYGBj8/Pw+fPjwfxSMglEwCkgHDP8HJRgSFeHr168ZGBgyMjKWLVs2WhEOIJgzZ460tDQjI6ODg8P/UTAKRsEoIBEw/B+UYEhUhGvXrmVgYOjv7182WhEONPj9+7eDgwMDAwM3N/f/UTAKRsEoIAWMVoTkg4KCAkVFxdmzZy8brQgHB1izZg0jIyMjI+MQdf8oGAWjYEDAaEVIPtDT03N2dl60aNGyZcuWL1q8wypyEB66PX/+/P8kgvnz57+kxqHbAwIePnzIy8vLx8fX1dX1fxSMglEwCogAoxUhmeDNmzcMDAxpaWnL8IKBrQhHLMjIyGBgYDA2Nv4/CkbBKBgFhMBoRUgmWL58OQMDw8SJE5fhBaMV4UCBxsZGBgYGfX39/6NgFIyCUYAXjFaEZIK/f//GxsbOmjVrGV4wWhEOILh+/ToDA0NYWNj/UTAKRsEowA1GK0LywcyZM0crwkEO5s+fz8DAkJ2d/X8UjIJRMApwgNGKkHyAXBEuX7DoAKf52vquZahgYCvCK1eulJB+CFlJScnVq1f/DxewYcMGZmZmQUHBZ8+e/R8Fo2AUjAIMMFoRkg9QKsLR7RODGOzfv5+Dg4OLi6uuru7/KBgFo2AUoILRipB8MFoRDi3Q0tLCwMDAyMiorq7+4MGD/6NgFIyCUfD////RipASMFoRDkUQFBTEzs7OwMAgLS395MmT/6NgFIyCEQ9GK0LywWhFOHTBwoULubi4ODk5u0b33Y+CUTDiwWhFSD4YrQiHOmhtbWVhYdHX19+wYcP/UTAKRsFIBQz/ByUYrQipAkYrQoLgwYMH8fHxDAwMGhoaa9eu/T8KRsEoGHlgtCIkH4xWhMMG7Ny5U1BQkIGBQUxMLC4urr29fdeuXVevXt2yZcvr16+/ffv28uXLO3fu7Nq1a/v27QsXLoyMjKyvr+/r62ttbX05ZM9lHQWjYBRAwGhFSD5AqQgXLd5pFDLYDt1+9uzZzJkz/5MIZsyY8fz58/8jD9y8eXPixImenp4MeAEfHx83N7e4uDg7OzsDAwMDAwMHB4eSklJ9ff2SJUt+/fr1fxSMglEwpMBoRUg+mIlUEeICA1sRjgLywM+fP+/evXv79u2bN2/u3bt3165dx48fv3z58sOHD79///4fCbx69Wr79u0dHR2GhoaCgoIMDAwMDAxiYmJGRkYyMjLMzMzS0tKSkpKurq4BAQExMTEZGRmurq6FhYU7dux4gGMLx/v373ft2lVbWxsVFcXExMTIyMjCwsLPzy8hIcHPzy8qKmpkZJSenj579uw3b978HwWjYBRQDBj+D0owWhGOgqEI7t27t3nz5hkzZtTW1iYlJWVlZdna2paXl2dnZyckJISEhHh4eEhISIiJiTEwMDAwMDAzMzMyMioqKnJxcbGzs7OxsbGysjIwMDAwMHBwcLi5uZWXl0+ePHnatGm9vb2tra21tbWBgYHq6uosLCwMDAwMDAysrKzW1tatra1RUVGHDx/+9u3b/1EwCkYBiWC0IiQfjFaEo4Bs8PHjx2PHjk2aNKmuri42NtbT07OgoKCxsbGrq2v69OknT558/Pjx48ePnyGBpzDw+PHjhw8fnjx5cvHixZGRkc7OzrKysgwwYGhomJ6ePnfu3MtUnej9+PHj169ff/78+X8UjIJhB0YrQvLBaEU4CigBv3//fvDgwVMYePLkybNnz54+ffqMFPD06dMnT548fvz4/v37R44cmTlzZlZWlrW1NQ8PDwMDAx8fn7i4eFpaWm9v77179/4TAs+fP9+9e3d5eXljY6OhoaGYmBgPDw8TExMDEuDh4WFjY+Pn5w8ICAgICHBzc5s1a9arV6/+j4JRMGTBaEVIPpiJVBEuX7R4v6DNmqbeZahgYCvCGzduVFdX/ycRVFdX37hx4/8ooDGAV4TPqAqePn36+PHjR48e7dy5s7OzMyoqSgw2EisgIKCpqWlqapqbm1tVVdXe3l5YWBgcHJyammptbS0sLMzAwMDAwMDCwqKhoWFnZxcZGVlcXFxfXz979uylS5cuXrx44cKFEyZMyMzMTEpKCgkJUVVVZWNjY2BgYGBgEBQULC0tXbhwIdV7jR8/fty+ffuqVasmT55cUVERFxcnKirq6uoaGhqalJQUHR3t7u4eGBhoZGTk4eHR0NAQGBjo5uZmYWFRWloaEhKSnZ09Z86cdevWLV++/O7du/9HwSjAAKMVIfkApSIc3T4xCkgENKoI0cDTp08fP358+/btTZs2dXZ22tvbq6mpWVtbGxgYqKqqcnFxsbCwODg4JCcnt7e3r1q16vz5848fP378+PGTJ0+ePsXitqdPnz59+hTCePLkyePHj69evbpkyZKwsDBdXV0GBgYGBgZVVVU7O7v29vYlS5YcPHjw3r17v3///k8E+PXr1/Xr17ds2RIeHp6Zmenk5CQtLc3AwMDAwMDIyCgpKWloaOjh4REUFGRiYmJiYqKrqysrK8vDw8PHx8fPzy8pKSkjI8PPzy8gIKChoaGsrMzLy8vHx8eD1K9lZmZWUVHx8fGJiYlJT0/ftm3bjRs3fo2u9R3ZYLQiJB+MVoSjgBLwmy4VITJ4+vTpkydPnjx58vjxYwgJAU+fPn1KsTOePn366NGjkydPTp061c7OTklJSUVFhZubmwEGuLi4REVFlZWVlZSUFBQU5OXl5eTkZGRkpKWlpaSkJCQkGJCAlJSUra1tXFxcQ0PDokWLjhw58ujRo8ePHz9+/PjJkydHjhyZMGHCggUL5s+fP3fu3Dlz5syZM2f27NlzYGDu3Lnz5s2bO3funDlzZs2aNXPmzBkzZvT19ZWXlyckJPj4+Jibm4uJibHBurMMDAyMjIzS0tKenp5xcXHFxcWdnZ3z5s3bvHnzyZMn79279/nz5/+jYPiC0YqQfDBz5mhFOArIB/SvCOkDnj59+vTp0ydPnjx+/PjRo0cXL17ctm3b3Llzm5ubk5KSqqura2tr6+rq6uvrGxoaGhsbm5qaWlpaWltba2trZ8+evXv37tu3bz9+/PjJkydPnz59ii18IBXh0qVLl5EFFi9evGDBAkgdOXHixOrq6tTU1JCQECsrKysrKz09PWVlZQkJCeRanIGBgY2NTVRUVFNTU0BAwMvLKzg4OCYmJiUlJTc3t6ysrL6+3tvbu7u7u6+vb+LEiVOmTJk2bdqECRMKCwtnzZo1C1YZT58+PTc3t6ura86cOTU1NRMnTlyyZMmyZctWrFixatWqNWvWrFu3bsOGDZs2berp6Tlw4MDp06evXbv28OHDN2/eoG3d+fXr1+fPn1+9evXgwYOrV6+ePn364MGD27dvX7NmzaJFi2bMmNHX19fc3FxZWZmSkmJsbJyfn19WVubl5VVbW5uRkbFo0aJVq1Zt2rRp/vz5x44du3jx4u3bt589e/bhw4ffxHXfhxMYrQjJBzNnzpw5c+bSpUuXLl26bO78g6wma2o6lqICtIrw0aNHT0crwlHw////4VsR4gFPnz598uTJs2fPnj179vz582cw8PTpUzhJDIBXhFQES5YsWbhw4YIFC+bNmzdnzpzZs2fPmDFj6tSp7e3tVVVV2dnZ8fHxISEhXl5e9vb2VlZWpqamBgYGOjo66urqSkpKsrKykpKS4uLioqKiwsLCgoKCAgIC3NzczMzMkOFZPj4+fn5+fn5+VlZWbm5uXl5eZmZmdnZ2Dg4OdnbozhkWFhYmJiZGRkYGHAAiy8bGxohXDScnJz8/v6ioqJSUlIKCgqKiopaWlpqampKSkpycnKioKDMzMycnJwtsEw4mYGNj4+Pjk5CQUFZW1tHRsbS0dHd3DwoKiomJSU9PLyoqqqmpaW9vnzhx4pw5c/Lz85djgNLS0mXLlkECdvHixYsWLVqwYMH8+fPnYfTUp0+fPm3atKlTp06ePHnSpEkTJ06cMGFCf39/b29vT09Pd3d3V1dXZ2dnR0dHW1tba2trS0tLc3NzU1PTf+qB0YqQfDBz5swJMDCps/sgq8ncnIoJGGD9+vUPkMDT0YpwFPz//x9WEY4CMsDBgwcnDHfQ29vb3t7e2NhYXV1dWlqal5eXkZGRlJQUExMTHx+fnJyclpaWlZWVn59fWlpaXV3d2NjY3t7e09MzgRTQ39/f1dXV0tJSV1dXUVFRVFSUk5OTlpaWkJAQFRUVEhLi5+fn4eHh5ORkY2NjZmZmYGCgpaWlqqoqLy8vKSkpLCzMy8vLhg0wMjKys7NzcHBwcHBwcnJycXFxcXFxc3Pz8PDw8vLy8vLywVoGAgICgoKCgoKCQkJCwsLCwsLCIiIioqKioqKiYmJi4uLiEhISEhISkpKSkpKSUlJS0tLS0tLSCgoK//79+08lMFoRkg8gYxFQsH3XQVaTgxPnHsQAt2/ffo8Evn79+p9e4NJoRTiIwe/fv9+PArLAkydPDo6CkQ2uXbv2n3pgtCKkDvjz6etBVpP3+0//H0xgtCIcBaNgFIwCgmC0IqQO+Pvz1wWHlG93Hv0fTODx48eTJk36TyKYNGnSk9Gr20fBKBgFIwaMVoSjYBSMglEwCkY0GK0IR8EoGAWjYBSMaDBaEY6CUTAKRsEoGNFgtCIcBaNgFIyCUTCiwWhFSB3w9+evk8q+n89d/z+YwO3bt8nYdtrY2Hjnzp3/o2AUjIJRMDLAaEVIHTC6fWIUjIJRMAqGKBitCKkDRivCUTAKRsEoGKJgtCKkDhitCEfBKBgFo2CIgtGKkDpgtCIcBaNgFIyCIQpGK0LqgNGKcBSMglEwCoYoGK0IqQNGK8JRMApGwSgYomC0IqQOGK0IR8EoGAWjYIiC0YqQOuDvz19nTaK+Xrv3fzCBBw8edHd3/ycRdHd3P3jw4P8oGAWjYBSMDMDwf1CCIVcRjoJRMApGwSgYomC0IhwFo2AUjIJRMKLBaEU4CkbBKBgFo2BEg9GKcBSMglEwCkbBiAajFeEoGAWjYBSMghENRivCUTAKRsEoGAUjGoxWhKNgFIyCUTAKRjQYrQhHwSgYBaNgFIxoMEgrQgUFhU+jYBSMglEwCkYB6eDnz5//SQGDtCLMysrSwwHU1dVlZWX1Bh8QExPT1dXVG0xAR0dHQkJCj0QgLi6uo6OjNwoGE9DV1RUTE9MbBaNgFBABSD0kcpBWhHjAunXrioqK/g8+oK2t/fXr1/+DCTx79szKyuo/icDCwuLFixf/R8FgAp8+fdLT0/s/CkbBKKABGK0IqQZGK8JRQDswWhGOglFAOzBaEVINjFaEo4B2YLQiHAWjgHZgtCKkGhitCEcB7cBoRTgKRgHtwGhFSDUwWhGOAtqB0YpwFIwC2oGhVxFev3791KlT/wcfWLZs2a9fv/4PJvD58+d169b9JxGsXbv2y5cv/0fBYAI/f/5csWLF/1EwCkYBDcDQqwhHAT7w9tzpm4OrMh4F5ILv754/f/781afR+BwFo4DGYLBXhN/v7Z/dUl5e3jJ7/63BVCL8+nTrxM6dO3eeGFSu+n+2VZaB3XE28XXhr5eHJgSaSUtLmwVOOPb2/ygYJODtsVYLXgYGBgYGBnblqBXER+goGAWjgHQwqCvCXzdnO/Iq+2SVl2f5KLMzsBuW7no0GEqEt/tKlXlF7eOzAg1F2Rl43SeeGyxVyIsVAQwMDETXhW/3FcjyWlQu3blzaWXuPKK0jALag183ZzvyOnaffvnr/69Pl6Y7sjN4LnnyfxSMglFAKzCYK8InSzwZwtdB65hfL093O7Iz8LoPeHn9ZEUAu/V0qCt+PdqcLMvAIFuwD+rOAQYnG0Q922cnE9cvvDpRV7T6GFzVr0+vnj9/9/3/KBhY8HFbMrt698X/UPDrWLUoQ8CKF/9HwSgYBTQCg7ki3FfAgFQe/P///9fNee68A10XvlgRwGA2/eZ/OPi4LZmdYbDUhR83x7KXHoRUz+yOs2/++nVz1/arOCq3fQUM1rPv/P////+vl7sqtXkZGBgYGBh4LVpHB0kHEtyf54iSwl6sCIDF0ygYBaOAFmAwV4QvVgQw6HZfRKn0ft2c7cjO7jh74KrCj5tjGdgL9iHZf3O6Wdb02Y7solm7Br4Cudit7jjv/n9oV5VdWVmW3bD15Mf/2MDFbnVIl/vJigBew9Lt158/f359e6khO4Nsw0kkD44C+oKrE3UZlLsu/IeCX7uyRivCUTAKaAkY/g9i8HFfgShmT+vt5lh29uRt2Mt2OoCP25LZkVz1dl+BrPXsO//vzLZmkG09+3+AwcfNseylB///////19VJZgwMDAzwYVxMcGe2NYNswb6L8xzV28/BFb3dlTU6Fjeg4Neje0iT4XdmWyOi4/vVJa2ja2dGwSigLhjUFeH//2/3FcgyYHQAD5ayQwt7OgG0tatvjzUYsjPwSrvFx7tJ8zLIQirFj5tj6dpwf3tudoI2LwMDAwOvdmArfBnRxW5169l3ft1cESbLbtiwbmkycgCieQTawZaVFbVGcfjJBlGGgn3/RwGdwK+bK7BFJRTcmW0NjY3vV6c7srM7dp97+38UjIJRQD0wyCvC//9/3VwRJsuAkvl/7StgF204+Z9O4O2+Allei9zZ65dOiLfgZWA3bDj08tevl1fWTYDUKPegM3B0HcH6dXO2I69j9/Gn3399enR8ojsvA7thw7G3/6H1cViYLIz/69HmzKjuc29xeOT/230Fsgyo/caPm2Pp55NR8PFgqaxs8rrrz58/v7690oKXgV25dNdLeGz831fAwFCwb7QWHAWjgGZg0FeE//////VyV6kyOwOvReW6C/fvX1iaLMvuOI9exfTHbcnsiHHDXy93FcgyMMgWIAqqX58+ff////+vm7Md2elWe6C3BX5d7NZlYGA3bDj29v+TFQHsEBYKwOMRSA+X3bB06Ynr9y+sq7TghXZyRwEdwIsVAYiRz///obEBiZn///9DKsLUSaO14CgYBTQDQ6Ei/P////9fL0/PTjAUZWdgFzVMWIJrGSQNwMFSdoaCff8R4MkST3Z2dgZoVfH1RL0yA7uosjQvg2wY/eZu7s9zRFtRe3O6mbWjIzt7wIon/9/eu/f2PzrA7xGU8J19+iW9PDIK/r9YEYDe/357rtuRHb5i6f48RwYGhtFacBSMApqBoVIRDhi4M9uaAdGR+v///8fNsdbTd812ZJctPfjx////kDNm6HTCzK+XT1/++v8fsrAQUoVBwNvNsbINJ99ujmVnj90MF0UGBD0yCgYIPFniiRKX/////w8ZsZaFVIVPVoS5j9aCo2AU0A6MVoSEwJMVAewMsIGqXy93FciKNpyEjEWiteNpDSDTpRYzb/3/Dy0oIaOZ108sLTVkhxSaVyfqooyz/f//6/v3X////x9MHhkFqODXyQZZBsSKJij4uK9AFHWfzigYBaOANmC0IsQAv16enp1gpq6uru5Wuf3e9///356b6M7LwMAuKi3KzsCuXLrv7f/////fnG6GOtRIWwCpBVE6DvBlo+yi9pWwpYb7CuDV2q9Hm/O0eRkYGBh4tVuOfv3/f1B4ZBRgAb9uznZkZ5BN3gyNxf///////3ZdODwuR8EoGAU0BKMVIRp4u69Alhdy/GarmywDe8CKJ////////emFgztRRkAvdqvTbQ8HlloQO/h1sVsXusfy48FSWdmw2cfvP79/fN0ueGdjYD0yCnCC71dnu/My8LpPhE3Qvt2VJYp8At4oGAWjgFZgtCJEAR+3JbMjDrP5daxaFK1N/uvRlQtPv0OOQkYopDF4vSaMgcFt4aP////////93vZKN2VRdl5pM/iqlrf7KgNzZ69f2uoGX1D7YkUAQ+zmj/8hADKRefDCI3j1NyAeGQV4wa9HuyoteBnYRe3jy+PtRXlhffZRMApGAW3BaEWIAvYVMKAME6LzH60K5GVgYGBg4NXOQxnGoi34dXO6NYNs6cGP/5+sCGBnN8yYvX7n0lYfZXZIL/HXy9NLct3U1d1yZx+GOerObGuGrF2//v/////JtmRZBihgh25RGyCPjAKC4Pu9E0snlJdPWHri3vf/o2AUjAJ6gNGKEAXcmW2NcrzpndnWKBXh////f316NQCXpf66OduRXba0tUAUaef7230Fsgw4Rs9+7StgFy3Y9/H/x23J7Nbtp1/++v//F+QGD6gHB8Yjo2AUjIJRMNjAaEWICt7ua61E6iGhbr57+/QpHRvpkAvK38Fs/HVztiM7A9pI7cfNsWj7CeHg18VuXQbZgn3rCmA9w////////3ZzLPvoWsRBCd4eax29FHIUjIKBACO+Ivz1aFdroJm6urpZYOu6Ky/RiqF9BfCa5+2+UmV20cL9n//TAbzdV6rMzsDAwMDAwC5q33ro5a//////ujnbz2/axY//EeDXrixYRYjpkV83Zzuys7Ozo1SEo2sRByv4eLLVkBd9D8UoGAWjgA5ghFeEb/cVyMqGTVy/c+f6iVHK7Axod/G9WBEAqTTe7itVZqfbRYhPVgTwOk6/+v3////f70FOn0Q5VuT7p0+//v+H1nS6E6/+//8fh0fe7itVZmdghxr2/////7+OVYuKNiAOZxsFgwh8/w6LplEwCkYBHcHIrgivTtRFGib8BTnTVLZg39v/EHBntjWD9exz+0rpWAtCjtRCLPj8//8X5GpBWF14f7E7O7uosrq6sig77Fg3PB6BbB5kV/apnL1+aauPMjtEeBSMglEwCkYBBIzsivBgKTta7wh64nHpwY///////3FzLAM7Ly89a8H/EFuh450w8OvlrgJZBujpad+fXliHtq6QgEe+39s+Id5NXd0tfsJ2uJ5RMApGwSgYBf//j/SK8OO2ZMx75d/uK5CF1SBPlngy0LcW/P8f4iyME7d+3ZztiPNC4o/bCHhkFIyCUTAKRgEOMLIrQsjSSsw652SDLPTETtS7wmkHvkMOfDl44en3/////78zz5GdAX0M89e5dnWUMVMkQNAjo2AUjIJRMAqwg5FWEb49tyTXTV1dXd0tF3Kl7q+bsx3ZGWBnUcMAfa9of7uvQBY268fArpy57t73/2/3FcgyoF0riLTg8xfGneaDwSOjYBSMglEwBMGIqgh/3ZztyGtYunTnzp1LKy14GdiVo1bc/PUfevtbGOKaw7frwhkgqzHpAD5uS2aH7ZL/9fL07DBZaAUIneYLgx6k9uvmbNga0Y8HS7HdaT7AHhkFo2AUjIIhCUZSRfhxWzLyHfKQU44hw4m/bq6AbDrInb1+/ewMQ3ourdxXwKA/Camugi6MMWw49vb/d8juCQZeaXVpXgbYGtEXK3DeaT6QHhkFo2AUjIIhCUZSRXh/niNK/fH//6+bK8JkGaDnVP96eXpJZaCZullg5ZJzdKw87sy2ZvBc8uQ/Mni7r0CWAbb65fu76yd27jxx/R2sn/diRQADUo3+//9/WK+24eSvgfPIKBgFo2AUDEkwkirC/1cn6jLABiFh4NfN2Y7s0KpwgMCTFQHskI4pEvh1c7o1rnNE/z9ZQuBO81EwCkbBKBgFxIIRVRH+f7s5lh2zBnmyIgB9Fx59wa+bsx0x3fVrVxYDLmf9Otkgy4BReY7eaT4KRsEoGAWkg5FVEf6HdJuQjxz7/x/rUCO9wdt9pcrskCUv/2HgZIMo7htzoZXn6J3mo2AUjIJRQCEYaRXh///QA8sMGw7B6pxf59rVce1UpyOATPPxWlRuv/XpF+TGXPwDnZDVPqN3mo+CUTAKRgFFYORVhP//////9lirBS8Du7JP69KdkGvd0QYZBwj8enl6dpQy5NoJduVMlN4eVvBr9E7zUTAKRsEooBAM/4rw7cuXWKuT74jzN3cRrHCoDX7dXDEbZ5WF48ZcXB75P3qn+SgYBaNgFFAAhnlF+PZctyO7ct1xulwiSCz4dXOeOy87aZ23QemRUTAKRsEoGA5gOFeEb891O7KjL4z5/////1+Pdq07RkItRE2AWgt+f/r07X8oeLuvVNsd07H//w9Oj4yCUTAKRsEwAcO2IkStPN7uK7WPWgNdTvl6SwIvg3L9ia//6Q3QasGr0x3Z2aM2vP7////////fnut2d8ecq3w7GD0yCkbBKBgFwwcMz4oQo/JAu1iXXndKoABstSDssl1cYFB6ZBSMglEwCoYVGI4V4f3F7vgqjwECX49WKpNWCw5Sj4yCUTAKRsHwAkO/Ivz16dGFHbNbyssnLIUdx/n22oV7EMZgrDyIrAX///8/yD0yCkbBKBgFwwIM7Yrw16PNmcrsDPC7/BgYeC1aYRvlB0Xl8evl6SW5bupmkDsD/////3pDFOFa8NenV8/fff8PBYPCI6NgFIyCUTBcwVCuCN/uK5BlNyzdfg9aZfx6eXp2mCwD7Lait+cmVg5w5fF2X4Esu7JPVnm8BS8DrCp7e+8evloQstmfgYGBHXbg2iDwyCgYBaNgFAxfMHQrwo/bktlhFxUhgbf7CmQZGHS7Lw6CiuP+PEeYCyGXDMLqQtzgyYoAdtnkddef3z8+0Z2XQbYAUheOglEwCkbBKKAVGLIV4cfNsTjuZvh1c7o1w4CfHXpnnrt7cqx6wb7/MECwLvz16N6jXQWIw0J/3Zw3WheOglEwCkYBrcGQrQjvzLbGedHCx82xDAyQKujXrwGqRH7dXBEmy4B2PyDeuvD1hih2Xl5es+k3/8PBaF04CkbBKBgFtAZDtiL8uDmWgR3b3Xu/Pr26PtebgUHU0cdNmpc9c+eAbTeH1IXQGUsY+PVy12wUASTw/ep0R3b4xfRQ8OvmPHdedp9lj/6PglEwCkbBKKABGFIV4a9HmyfCK5EXa8LYMe6b/3+hS5mBV1SUnUE0oHv9wQuPPqFK0xj8enl6doKZurpbJXSNKNa6EA18v3di586DV17++v//P7QuRLuj99fLl8jcUTAKRsEoGAVUBEOqIryz2J0XMbD46+Z0awYG9GHDty+fnuvWZdCdePU/ncGvR5uTMdeI/sdbF/66Oc+dl4GBgYGBAb6nAmtdOApGwSgYBaOARmBIVYSwigNWyfx6uatAloEBepft////f326tb3UkH1Arhe8M8+RF2bvy1WBDAxwZ/7/dXNFmCziUFE4eLs5ll02ed297/9/vdyeJTpaF46CUTAKRsEAgCFWEf7//x+1Lvz/9hz8KlsGBgYGBnblqCVXoRsLaQreHpvQinRz7sVudYasXb/+/////+2+Alnd9nVLwmQRzvz1/fuv/2jgZAP8PvlfF7t1ZS0sZJHrwtm5raM3S4yCUTAKRgHNwdCrCP//R68L/////u7+hYM7dx68cP8dPerA//////91c7Yjr+M8+LLVO7OtGUym3fj//848R3bIWtG368IZGBjYlRM23P+PBbxYEcCg3n3x////b/cVyIpm7Xr7/8mKAHYG5DMCRsEoGAWjYBTQGgzuivDXo815boFZLbPXH7xw/9WnX//hAKMupCN4e+3Ko1//////9evX/183N8yGLIyBiD5ZEcAO283/a1eWetOW/YdxOfHJmty8zY9+vd2VJQobBb0z21rWwt4+c/MjHHpGwSgYBaNgFFAZDO6K8M5id14GJMAuqmwGrxffXR2QuvD1jjRRBtlkaFV1Z54jO4L3//++Agb+6mP//////+vmbEdZ6MgnMvj16dXz50gnid6f58jgNP/B//////9/uzlWFushAaNgFIyCUTAKaAUY/g9uAOn5yYatuPTq+fUTO5dOKI93U5eG1o7svLzsDAwMvO6L4SOU9ABQN0Frv183ZyPVhVcn6jKwG0aWl8db8MK6eQgAP0eUgYFd1G/iubf/IRWhaNbm+8/vH+9GGWwdBaNgFIyCUUAPMNgrwv///0PrHdT9B9/fwetFswQUGboAqJugtR9yXfjr0a5KN3V1s4TZp1/++o8Mfl3s1mV37D798tf//9/vLQvnZ4AsJH27rxSy2ofXYnR1zCgYBaNgFNAdDIGK8P///9B6J4z+NR4W8OvloQmBZtBrnyC13///yHUhDLx9+RLB+f//YCk7ZG3M//9v9xXIsjtOv/r9+9Xt++78hwyXIkZLR8EoGAWjYBTQDwyNivD///+DpS78dbJBFtqv+/70+OwwWXjtB60LM3e++P////+3+0qVUTY0vlgRwGA2/eb///Ba8P//G9NMIGKjYBSMglEwCgYKDJmK8P///4OiLjxYyo5y6cXbfQWySHXhrnXnIKObT5Z4MjAwMDAg6sKrE3UZZJMrk2G14P///3/tymKA9RJHwSgYBaNgFAwIGEoV4f///2F14Zo7/wcI7CtAP+v717FqUcQYKRz8OlYtqp5WEysLrwvf7iuQZWAQzNj57v/////////1cl0sfK/FKBgFo2AUjIKBAUOsIvz///+vmytaB7BP+HFbMsqI5////+/MtlZ3dI+CdfPg4M5sa/bkzdc2J8Prwl83V4TJMjDwagdmlWcFavPSfe/HKBgFo2AUjAJ0MPQqwgEAvz7dOrFz586DF55+/w/t2MmGwQ5y+3VztiMv9muAP+4rENWdePXXI6S68P/3e/tnt5SXl5dPWAe9cWIUjIJRMApGwQCC0YqQAPj1clepMjuvtLq6sig7A7thw6GXv37dXBGlzM7AK62urq4syg6ZtXx7bMns/fcglSMcXJ2oy5687eN/1LpwFIyCUTAKRsHgAaMVIX5wZ7Y1fG3L9+tT7BnYoVvef326dWL97Amz15+49/3//48HK5XZGRgYGNiVoyYiV4cftyWzW8++8///aF04CkbBKBgFgxOMVoSo4O2541cR1dj/m9PNoHdK/Hq5q0CW3XH2zV9vz01sXfPkPyqA7JYQDaxrT9DmZWBX9mmFHZx9daIu/I6JR5uTlR1nj9aEo2AUjIJRMIjAaEWIBD4eLJVlYID3AP9DLleK3fwRUQv+///rWLUoQ+xmzDlBSF2oXLrv6dPjS3It4NXhx23J7J5LoDXnr1+//o+CUTAKRsEoGERgtCJEBlcn6jKws7Mj6sJfx6pFGWQtLGC14P//kMoR0kvEAPC68O3//9/h1WFYrJsgrFM4CkbBKBgFo2CQgdGKEAV83JbMbtY0vUAWXhe+3Vcgy8Bg1nT85a////////VoXayobMPJX/9xAOS68P///79eXllXacHLwO6++P7/UTAKRsEoGAWDD4xWhKjg18VuXdHqQ093FcjC6sJfjzZnKrMzMPBKq6tL8zJA1ojiA2h14f///399+vT9/ygYBaNgFIyCwQhGK0J08HZzLLvnkie/XiLVhf9/fbp1YudO+FZCggCzLhwFo2AUjIJRMDjByK4I3+6rDMxdchy1cvt1skFWtPrYr/+odSGp4O2+UmV23oQtr/+PglEwCkbBKBjEYCRXhBe6lBkYGNjZ2Rl4LVCqwycrAngDVjz5/5/SuvDpU3K0jYJRMApGwSigIxjJFeH/t8caDNkZZFPnLmv1UWZn4NVOmA2pDn+dbJCFnYb96+Wu0qjZBKYFR8EoGAWjYBQMVTCiK8L/8LqwYN/b7/f2T4yCV4d3l3iyx24eneEbBaNgFIyCYQ9GekX4H7ku/P////d7+ycmaPMy8BoaqjLoTrz6fxSMglEwCkbB8AajFeH//2h14f///78/PT47QVvUcXREdBSMglEwCoY9GK0IoQC9LhwFo2AUjIJRMDLAaEWIAKN14SgYBaNgFIxAMFoRooC3xxoMeaEXLY2CUTAKRsEoGAlgtCJEB6P3Q4yCUTAKRsGIAqMV4SgYBaNgFIyCEQ1GK8JRMApGwSgYBSMajFaEo2AUjIJRMApGNBitCEfBKBgFo2AUjGgwWhGOglEwCkbBKBjRYLQiHAWjYBSMglEwosFoRTgKRsEoGAWjYESD0YpwFIyCUTAKRsGIBqMV4SgYBaNgFIyCEQ1GK8JRMApGwSgYBSMajFaEo4BicGe2NQMMWM9GOagVWQpVcl8BA0PBvv/YAW5t2MC+AgYEIKSYIICZVrAP7gzr2UsKGBgKcLn2P8TBlNi8bx92w+/MhrmAfLMJgX0F+P1GMbgzGxY4CNZgAsS4iuaBNAoGFIxWhKOAMnBntjVSEbGvAKnQRuH8h/DhSvcVINioYF8BHm1o4M5sa1RJfIqJAXdmIwrFfQUUGUU0wG0PbhkqAppbghymgxEQ4z6aB9IoGFAwWhGOAooAeiGCxMdSdiAJITFRARYJLEL///+HWIYhgVWQaIBk1Z3ZcK/QFiDZiQro4wKc1lML0Mcb5ANi3EfzQBoFAwpGK8JRQFWAVKjsK0Dt2aGCfQU4SpZ9Bfi0IQMku3CDfQUMMICqFlPizmxrBmygYN//fQUMDAUw1yKrQ9IKZf6H8GAARZf17H2z4VIQmX0FDDAAEYABJAmGgn3/9xUwWM+eXcDAwMDAUABRuK+AAQYQdu8rYGAomA23xXr2nf93UHgYAE0HzPT///8j6WRgYECW2VfAAAcIUVQNcPE7s62hNqOzMEIDAu4ghFHEEWBfAQOym61n4/IlQpQBxSgk8YLZs62R9CDJIOnYV4DMGwXDDYxWhKOAmuDObGtEgXFntjUDFMCEEAB3yYJXGzLAbQQU3JltjSgWkTl3Zlsjad1XAJf4v68ALnFntjUW4TuzrRGq78y2hojfmQ1Xu68AIf//zmyoAigTLnNntjVM4v++AjgTFdyZDTf1/74CBoRuiBSci8zZV8DAADcOwoHK/N9XgGAjgX0FDAxwHf/3FcAV7StACEN4EO6d2Qhn/d9XABWFMOHid2ZbwyTuzIapR2UhVN+ZbQ1TjNMUFLCvgIEBLgHhwPTsK4Cz9xUwwBWhcPYVwNX831fAwADn7SuAM1Et31cAZ46CYQhGK8JRQD2wr4ABo7jYV8CABOCFzP99BRhKkcG+AgYkgNCGDPYV4JCAgjuzrVEUwPl3Zluj2r2vAOYYBOv/ndlQ5f+RhJEFkQBCeB+ao+Ayd2YjitX/EK717Dv/////v68AWQIJIKn5v6+AAc5Glfn/H5mPqm5fAZLRCEUoAEXNf0w+DMDF9xUw4DAHRRRuHZyBxoKY9v//fwgXIoHTFBSwD8UJ+wqQzEJoQBH+j+AjVEDAvgKYYfsKoAwoQKjcVwDROwqGJxitCEcBlcC+AgYGaKmBHdyZbc3AAC9O9hXAmfgBqjZkQMgIDPk7syEV4L4CBkxQsO//f4gUhPX/zmxruH/gwnAGKoCrvTPbmgEDwGWsZ9/5DwVIXBxmoqhBU4TK+w9RWrDvP7oMCu/ObIRxSGBfAaoomrI7s5F8BDELSQgi8P8/iiACWM++A5Gxnn3nP3bW////kbh3ZuMyBQXsK0C2GoV3ZzaSUdaz7/xHgH0FDNaz7/zfh6z8P5LCO7NxWo6uZxQMLzBaEY4CaoB9BQywIgMfuDPbGlae7CuAsQgDJG3I4M5sa6x27iuACO9DtwIqcGc2VuP+//8PV/L/P6rxcGE4AxXA1e4rYIAwMABcCQQgcfcVYDUTRQ2aIlTefyQBOOP///9oPGTjkMC+AlRRuLJ9BQwMDAxwA/YVwJlQsK+AAQIK9kE41rPv/McEcAOxsv7//4/ExW0KCthXALEUAlB4cKPgDBiACaAo/48Q/7+vAKfl+wpQ9YyC4QVGK8JRQDHYV8CAWUrcmY2tttlXABNEsFAAIW2oAJs4woQ7s61RyjU4f18BMQUeXPl/JGFkQSQAF74zG2Y7OrgzG6oEApC4cMPRAZIaNEXIMv//I/NR1aHwEIpQAIqa/wg+jIaCO7Nx+AyqDqf8/zuzYfZiY/3//x+Je2c2LlNQANROKEDh3ZkNNQpV+D+Cj1ABAfsKoMnhzmyclu8rwCUzCoYDGK0IRwFlAHfhsa+AAVrAwMCd2YgSaF8BWdrQwJ3Z1gwo5uwrYEDw78y2RhiFzEFm/4dogvL2FcC135mNsBYhfGe2NYz5H8Kznn0HTv+HKEVVAOHdmQ1X8v8/CheJiQqQJfYVQM2BgDuzreFuRuGgqkPh3ZmNMA4J7CtgYECo2lcAMwrB+g/hMEBV7StAkvi/rwAq/H9fAQOM+R9iG4R3ZzbMXmys////o3BxmYIC9hUgi6Lw7szGbhQyZ18Bwgf7ChgY4Lx9BQxwRRCjoLx9BTDWKBiOYLQiHAUUgX0FDJgAUWKgSSMk0GUYGBiQpNEk4eJYwZ3Z1gwIAC3R4ADJKBQpFF0ImX0FcOvuzIaXqMjCEAkYgCq4MxuhFqIaDmDaUJWgcO/MtmZggJuFAHdmI9TsK0AYBQX7ChhgAKYKIohQh8JDNg4J7CtggGwhgAAkFXfgggwMBfvuzLZmYIAYdme2NQMcIGn4v6+AAQEgaiGqIYqwsf7//4/OxWoKCthXgCyBwkM16s5shEvhSv7/h+iBgILZs5F1ICQYGBBa9hUgcUbBsAOjFeEoGAWjYBSMghENRivCUTAKRsEoGAUjGoxWhKNgFIyCUTAKRjQYrQhHwSgYBaNgFIxoMFoRjoJRMApGwSgY0WC0IhwFo2AUjIJRMKLBaEU4CkbBKBgFo2BEg9GKcBSMglEwCkbBiAajFeEoGAWjYBSMghENRivCUTAKRsEoGAUjGoxWhKNgFIyCUTAKRjQYrQhHwSgYBaNgFIxoMFoRjoJRMApGwSgY0WC0IhwFo2AUjIJRMKLBaEU4CkbBKBgFo2BEg9GKcBSMglEwCkbBiAajFeEoGAWjYBSMghENRivCUTAKRsEoGAUjGoxWhKNgFIyCUTAKRjQYrQhHwSgYBaNgFIxoMFoRjoJRMApGwSgY0WC0IhwFo2AUjIJRMKLBaEU4CkbBKBgFo2BEg9GKcBSMglEwCkbBiAajFeEoGAWjYBSMghENRivCUTAKRsEoGAUjGoxWhKNgFIyCUTAKRjQYrQhHwSgYBaNgFIxoMFoRjoJRMApGwSgY0WC0IhwFo2AUjIJRMKLBaEU4CkbBKBgFo2BEg9GKcBSMglEwCkbBiAajFeEoGAWjYBSMghENRivCUTAKRsEoGAUjGoxWhKNgFIyCUTAKRjQYrQhHwSgYBaNgFIxoMFoRjoJRMApGwSgY0QAARVvgXw1PcJ8AAAAASUVORK5CYII=",
      "text/plain": [
       "This front-end cannot display the desired image type."
      ]
     },
     "metadata": {
      "image/png": {
       "height": 436,
       "width": 600
      }
     },
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "\n",
      ". global stata_kernel_graph_counter = $stata_kernel_graph_counter + 1        \n",
      "\n",
      ". `5a6d50033241a343dc684724381211ea'\n",
      "\n"
     ]
    }
   ],
   "source": [
    "qui graph export \"figure.png\", as(png) replace"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Stata",
   "language": "stata",
   "name": "stata"
  },
  "language_info": {
   "file_extension": ".do",
   "mimetype": "text/x-stata",
   "name": "stata"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}
